How to calculate the degree of spin polarization in ferromagnets. 
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Different ways to define and calculate the degree of spin polarization in a ferromagnet are dis- 
cussed, particularly with respect to spin-polarized tunneling and Andreev reflection at the boundary 
between superconductor and ferromagnet. As an example, the degree of spin polarization for differ- 
ent experiments in Fe and Ni is calculated in the framework of the local spin density approximation 
(LSDA) and used to illustrate the difference between various definition of spin polarization. 



Although solid state physicists use the notion of a de- 
gree of spin polarization (DSP) of a ferromagnet (FM) 
rather often, it is not well defined. While the total mag- 
netic moment is uniquely defined as the difference be- 
tween the number of spin up and spin down electrons, it 
tells us little about how much do electrons with the dif- 
ferent spins contribute to transport properties. In view 
of the growing number of experiments probing various 
aspects of spin polarization Q, it becomes increasingly 
more important to be able to calculate the DSP in the 
framework of the conventional band theory (and eventu- 
ally beyond it). Importantly, the DSP can be defined in 
several different ways. Therefore it is crucial, in order to 
compare the calculations with the experimental data, to 
make sure that a proper definition of the DSP is used. 
In particular, spin-polarized tunneling in various forms 
H, including Andreev reflection provides valuable 
information about the spin dependence of the electronic 
structure, but this information may be obscure and not 
very useful unless the measurements are backed by the 
calculation appropriate for the experiment in question. 

Let us consider an extreme example, the so called half- 
metallic magnets, for instance Cr02- Such systems do 
not have any electrons at the Fermi level in one of the 
two spin channels; they are metals in one spin channel 
and semiconductors in the other. Half metals have 100% 
spin polarization according to any sensible definition. On 
the other hand, for a regular magnetic metal, which has 
Fermi surfaces in both spin channels, it is not obvious a 
priori how to define the degree of spin polarization. 

The most natural definition, and probably the one 
most often used is P = (JVf - N l )/(N 1 + N L ), where 
A r |(j) is the density of electronic states (DOS) at the 
Fermi level, defined as (h = 1) 



Ni = 



(2 



^3 E / ^ m )0?k 



1 



(27T) 



W— , (I) 



and E(v)k a i is the energy (velocity) of an electron in the 
band a with spin i(=t or j.) and the wave vector k. A 
typical experiment that can probe P/v is spin-polarized 
photoemission. This definition of the DSP may be called 
"iV" -definition, P/v, and its usefulness is limited by the 
fact that the transport phenomena usually are not de- 
fined by the DOS alone. This is particularly true for 



materials which have both heavy d-electrons and light 
s-electrons at the Fermi level (a good example is Ni). 
While the DOS is mostly defined by the former, the elec- 
tric transport is primarily due to the fast s-electrons (cf. 
a semiempirical recipe of defining DSP via partial s-DOS 
in transition metals, ^)). 

Classical Bloch-Boltzmann theory of electric transport 
in metals Q lets us separate the current of the spin- 
up electrons and the current of the spin-down electrons, 
and to define DSP via the current densities Jfm a s 



(J T - J;)/(J t + JJ, J 



TU) 



oc 



(Nv 2 ) 



T(1) T TU)- 



Assuming 



the same relaxation time r for both spins, this definition 
leads to the "Nv 2 " DSP: 

JW = «^ 2 ) T - <^ 2 ) i )/«^ 2 ) T + (Nv 2 ) ^ (2) 

where (iVu 2 )^^ is defined as 

(Nv 2 ) i =(2^)- 3 ^ I v 2 m 5(E kai )d 3 k 

= (2tt)- 3 ^ f v km dS F . (3) 

a J 

This quantity is sometimes denoted as y^) s & n d is pro- 
portional to the contribution of the corresponding elec- 
trons to the plasma frequency (see, e.g. Ref. B). If 
spin-dependent or spin- flip scattering is present, the to- 
tal current in each spin channel depends on the charac- 
teristics of both spin subsystems, and the expression for 
the DSP becomes very complicated. 

Unfortunately, it is hardly possible to measure Jj 
and J| separately. A typical experiment involves spin- 
polarized tunneling between a FM and another mate- 
rial. In particular, it is possible to measure tunneling 
currents separately for both spin polarizations for a fer- 
romagnet/superconductor contact. The question arises, 
whether the DSP measured in such a way is P/v or P^v 2 ? 
To answer this, let us start from the simplest case, a 
ballistic contact with no barrier, and neglecting possible 
mismatch of the Fermi velocities at the contact. We re- 
peat the original Sharvin j^] semiclassical derivation, but 
allow for arbitrary Fermiology. Following Sharvin, we 
assume that an electron going through the contact ex- 
periences the acceleration by the electric field so that its 
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energy increases by eU. Assuming that the field changes 
the electron's quasimomentum from Kk to Kk' , we find 
that the phase space for this process is defined at T = 
by the factor 

0(£ k ,)0(-£ k ) = e(E^ + eU)6(-E k ) = eU6(E k ). (4) 

The fraction of electrons with a given k that can reach 
the contact in a unit time is v x A (the contact plane is 
perpendicular to x and A is the area of the contact). The 
total current is 

I = tz To - / v x 5(E ka )dk =e 2 UA (Nv x ) = 
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where S x is the area of the projection of the Fermi surface 
onto the interface plane. For a Fermi sphere this reduces 
to the Sharvin result. Correspondingly, we arrive at the 
third, ballistic definition of spin polarization, Pn v . 

The next simplest model is that of a specular (<5- 
function) tunneling barrier with an f/s Fermi velocity 
mismatch. Here we need to take into account, in addi- 
tion to the (v x /v) factor, a finite barrier transparency. It 
depends (nonpolinomially) on the Fermi velocities 
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where v s m is the Fermi velocity in the superconduc- 
tor(ferromagnet), and W is the strength of the barrier, 
V(x) — W5(x) (for a one-band isotropic material W is 
related with the dimensionless parameter Z of Ref. as 
Z = W/vf)- Tunneling current is thus proportional to 
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In the large W limit this reduces to (Nv 2 ) and the mea- 
sured DSP is P/v„2. However, in the high transparency 
limit one cannot give a simple answer. 

It was recently suggested |§ that Andreev reflection at 
the interface between a FM and a superconductor (SC) 
can be used for direct probing of DSP. The idea is simple: 
the Andreev reflection can be visualized as two currents 
of electrons with the opposite spins flowing inside the 
normal metal towards its interface with a SC. At the in- 
terface (more precisely, within the coherence length from 
the interface) the two currents recombine creating the 
current of Cooper pairs. One can also say that to form 
a Cooper pair an electron has to leave behind a hole, 
creating an additional hole current flowing in the oppo- 
site direction. In a paramagnetic (or antiferromagnetic 
with time reversal symmetry) metal both currents are 
the same, so one observes in the superconducting state 
the 100% increase in the net current over the normal 



state. De Jong and Beenacker || suggested that in a FM 
the total Andreev current is defined by that spin channel 
where the normal-state current is smaller, because the ex- 
cess electrons in the other channel will not find partners 
to form Cooper pairs with. This statement was quanti- 
fied in Ref. Q via the number of spin-up and spin-down 
conductance channels, which they denoted as N^^, thus 
arriving at an expression for the ratio of the current in 
the superconducting and the normal state as 



I s /I n =Amm(N h N l )/(N ] +N l ). 
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The number of conductance channels cannot be directly 
evaluated. Besides, this expression can mislead the 
reader into a belief that the DSP measured through An- 
dreev reflection is Pm (i.e., defined by the DOS). 

In fact, Andreev reflection at a FM/SC contact has 
been recently attracting substantial theoretical interest 
This interest so far concentrated on such aspects as 
the symmetry of the order parameter, Fermi velocity mis- 
match, and generalization of the BTK formula onto 
spin-polarized case In terms of electronic struc- 

ture, however, all the work was limited to the parabolic 
bands/spherical Fermi surface model. While revealing 
important fundamental physics, such approach is of lim- 
ited practical importance, because in real ferromagnetic 
material this approximation is absolutely inacceptable. 
In this Letter we, on the other hand, focus on the band 
structure effects in spin polarization, and in this context 
we need a better definition for the "number of conduc- 
tance channels" . 

Comparing Eq.|| with Eqs|5] and [?], we observe that 
the DSP for Andreev reflection should be defined as ei- 
ther P/v„2, for a large barrier and/or diffusive current, or 
Pnv, for low resistance ballistic contacts. This is, how- 
ever, only the first approximation, while full expressions 
should include Fermi surface averages of more compli- 
cated functions of vp. For one particular case, a fully 
ballistic (Sharvin) Andreev reflection, an explicit for- 
mula, reflecting the physics suggested by De Jong and 
Beenacker, can be derived, which is both suitable for 
band structure calculations and also quite illustrative. 
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FIG. 1. Fe bands in the two spin channels. The 
linewidth is proportional to the partial s-character in each 
state. 
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FIG. 2. The same as Fig|l[ for Ni. 

De Jong and Beenaker treated incoming electrons and 
reflected holes as two separate currents. In purely ballis- 
tic regime, however, one has to take into account energy 
and momentum conservation (parallel to interface) for 
each individual reflected hole, so that Pjvu, with its inde- 



pendent averaging over each spin channel, does not nec- 
essarily correctly describe observable polarization. This 
statement may be quantified as follows: Consider an in- 
coming electron with a momentum k = (kii , k x ), which is 
reflected as a hole with the momentum q = (qii , q x ) in the 
other spin subband; qn = kn, = Eq±. This fixes for 
each k a countable number (one, two, etc) of wavevectors 
q satisfying this condition. For simplicity we will assume 
for the moment that there is only one such q and will 
denote it q. This adds an additional constraint to EqJ|, 
so that instead we have, omitting summation over the 
band indices a, 
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This can be rewritten in a symmetric way like 
I = (2ir)- 6 e 2 UA f d 3 kd 3 q f d 2 a[5(q|| - a)£(k|| - a)]M(£ q )]M(£ k )], 
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where we introduced a 2D vector a, transformed the 
(^-function of q into a ^-function of _E q , and introduced 
u =dEq/ (9q. We can further rewrite Eq|)| as 

I = e 2 UA J rf 2 aF T (a)F x (a) ={2-k)~ 3 e 2 U AS , (10) 

where F(a) = (2tt)- 3 / „ 5(k\\ - a)v x S(E k ) is half (be- 
cause v x > 0) the number of crossings of the line ky — a 
with the Fermi surface for a given spin. Using again 
projections of the Fermi surface for either spin onto the 
interface plane, one can define S as the overlap area of 
the spin-up and spin-down projections. This current can 
be expressed in terms of a spin polarization as in EqJ|, 
thus giving yet another, "ballistic Andreev" definitinion 
of spin polarization, similar, but not the same as the 
u Nv" definition, Pb.A. ^ Pnv (they are equal only if the 
Fermi surface projection for one spin is entirely contained 
in that for the other spin). 

So we observe that an experiment would probe differ- 
ent DSP's depending on the length scale of the prob- 
lem, which is defined by the size of the contact and the 
length at which the voltage drops, and how it compares 
with the mean free path. Furthermore, the transparency 
of the barrier can also influence the measured DSP. In 
the pure ballistic limit the DSP is related to the aver- 
age Fermi velocity, while in the purely diffusive regime 
it is defined by the average squared Fermi velocity p0[ . 
One may ask, why DOS is so often used as a measure 
of spin polarization, even though such definition is irrel- 
evant for transport properties? The answer is that the 
most common way to perform tunneling or similar ex- 
periments is to follow the details of the contact conduc- 
tance as a function of applied voltage. Probably the most 
spectacular and fruitful application of this technique is 



the tunneling spectroscopy of superconductors. In such 
case the characteristic scale for the voltage change is the 
superconducting gap and Debye temperature. The nor- 
mal state electronic structure does not change over such 
a small energy range, so the only important factor is the 
variation of the superconducting DOS with energy. The 
normal state DOS and velocity can be assumed constant 
and factored out. Of course, it is not the case when two 
different sheets of the Fermi surface, as in Ref. ||, or two 
different spin channels, are compared. 

Importantly, P^v 2 > Pnv and Pn are entirely different 
in real materials (Figs. ^,||). The reason is (and the 
"s-DOS recipe" works for the same reason) that in tran- 
sition metals one can often distinguish the pieces of the 
Fermi surface that are predominantly d in character and 
the pieces that are mostly s. The former have low Fermi 
velocity and, according to Eq.|l|, are responsible for the 
most of the DOS. The latter have high Fermi velocity 
and provide the main contribution to Eq.||. The larger 
the anisotropy of the Fermi velocity (angular anisotropy, 
in principle, works in the same way as interband one), 
the larger is the difference between Pnv 2 i Pnv and P/v . 

To illustrate this difference in real materials we present 
here Linear Muffin-Tin Orbitals JllJ] LSDA calculations 
of the corresponding quantities in Fe and Ni. It is instruc- 
tive to start from the band structure itself (FigsM^,||). 
The "fat" bands in these figures correspond to the states 
with substantial sp character. One immediately notices 
a qualitative difference between Fe and Ni: In the former 
all bands in both spin channels are heavily hybridized 
at the Fermi level, and one cannot convincingly classify 
bands as "predominantly sp" and "predominantly d" . In 
the latter the d band is so deep that one can single out an 
sp-like pocket in the spin-up channel (in Fig. 0, the band 
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crossing Fermi level between T-H) , and e?-like pockets (in 
Fig. near H). The Fermi surface in the spin-down 
channel is entirely sp-like. This is similar to paramag- 
netic Pd jl2j , which is the only Ad metal where transport 
properties can be described by the s — d scattering model. 
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FIG. 3. Density of states, N(E) (solid line) and 
(m)eff ( S ) = ( Nv f)(E) (dashed line), for Fe (left panel) 
and Ni (right panel). 

In general a large Pn may be due to any of the two rea- 
sons: either the areas of the up- and down- Fermi surfaces 
are different, while velocities may be similar, or the areas 
are not too different, but the Fermi velocity for one spin 
channel is much smaller than for the other. In the for- 
mer case the additional factors of vf or v F may change 
the DSP somewhat, but qualitative changes, or, in an 
extreme case, the sign change, are unlikely. If, however, 
"light" and "heavy" electrons are present, Pn is domi- 
nated by the "heavy" pockets, and Pnv 2 by the "light" 
ones, so the two DSP's are likely to be very different and 
possibly have opposite signs. 

The first situation is realized in Fe. One indeed can 
see that N(E) and (Nv%)(E) behave similarly (FigJ|). 
Correspondingly, the difference between different DSP's 
is only moderate (Fig. ||). On the other hand, in Ni 
the DSP essentially drops to zero when the factor v% 
is included. Interestingly, most experimental results 
indicate that the DSP observed in tunneling is positive 
and not too small (> 20%), in other words, the effect 
described above appears to be even stronger in reality 
than in band structure calculations. This is also to be 
expected: the local density approximation (LDA) has a 
tendency to underlocalize d electrons in transition met- 
als. For instance, in Cu the fully occupied d-band ap- 
pears in the calculations about 0.4 eV higher than in 
experiments, and is also too wide Similarly, LDA 

underlocalization of the d electrons in Ni leads to an over- 
estimation of the d bandwidth and of the exchange split- 
ting (by approximately a factor of two). As a result, 
the separation of carriers into sp-like and d-like in Ni 
should be even more pronounced compared to LDA cal- 
culations, and thus the effect of Fermi velocity on DSP 



even stronger. This leads, in turn, to the DSP sign re- 
versal, observed in tunneling experiments. 
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aspects of the spin-dependent Andreev reflection. This 
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FIG. 4. Degree of spin polarization for Fe, calculated 
as Pjv, Pnv, and P Nv 2 
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FIG. 5. The same as Fig§, for Ni. 
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